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Abstract—Investigation on time-domain responses of the 
monopole in the presence of an in _ intentional 
electromagnetic pulses (EMP) is carried out in this 
conference paper. An efficient hybrid numerical method 
based on FDTD method combined with the thin-wire 
model using subcellular extension is implanted for 
studying an EMP coupling with some typical monopole 
antennas widely used in the communication systems. Our 
developed program code is used to capture the responses 
of the monopole antenna. The transient response 
characteristics of the antenna illuminated by different 
incident EMPs are shown and compared in detail. These 
results are useful for further evaluating electromagnetic 
environment and electromagnetic protection of some 
complex communication. 
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I. 


In the past one century, the monopole antenna has 
been seen widely spread in the short-wave and ultrashort 
wave communication system for different platforms, 
especially for the close quarters communication and 
telecommunication between the ships and fleets. 


On the other hand, it should be understood that 
electromagnetic environment around us is now 
becoming more and more complex, because the number 
of electromagnetic interference (EMI) sources existed is 
being increased rapidly. Actually, most of 
communication systems are very sensitive and 
susceptible to EMI and in particular to an intentional 
electromagnetic interference (IEMI). 


INTRODUCTION 


IEMI refers to “intentional malicious generation of 
electromagnetic energy introducing noise or signals into 
electrical and electronic systems, thus disrupting, 
confusing, or damaging the systems for terrorist or 
criminal purposes” [1]. Such an intentional interfering 
electromagnetic energy can be generated by a non- 
explosive high-power microwave source [2], and even a 
very powerful ultra-wide band pulse antenna. 
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The IEMI usually works in the following way:(a) an 
intentional high-power or  ultra-wide band 
electromagnetic pulse is suddenly generated and it is 
incident on to antennas, metallic enclosure with slots, 
and transmission lines of the communication system 
through direct radiation or indirect capacitive or 
inductive coupling; and(b) as a result, it then causes high 
voltage, large current and rapid temperature rise or 
thermal stress in the inner circuit which may result in 
electrical or thermal breakdown and even destruction of 
the system[3]. 


With respect to IEMI coupling paths with most 
communication systems, we should at first take radiation 
coupling between the IEMI and the antenna into account. 


Further, we know that the analysis of antenna 
characteristics in the time domain has been carried out 
extensively in the past two decades for the sake of the 
time domain method can effectively capture wideband 
information and handle nonlinear problems. Among 
various numerical methods, time-domain _ finite 
difference (FDTD) is frequently implanted to solve 
transient radiation and scattering problem. Conventional 
thin-wire FDTD techniques can be classified into four 
approaches: the auxiliary differential [4], the contour- 
path integration [5], the weighted residual interpretation 
[6], and the modified of assigned material parameters 
[7]. These techniques have been improved by arbitrary 
orientation [8], cylindrical field behavior around the 
wires [9] and the concept of effective radius. 


Recently, a resistive thin-wire FDTD model with 
both the insulating and surrounding media being lossless 
and loss has been proposed [10], and the antenna excited 
by the coaxial lines is simulated by an equivalent model 
without any additional transmission line (TL) cell 
modeling. 


In this paper, we first improve the thin-wire model 
using subcellular extension [10] in the orthogonal 
coordinate system, and then a numerical procedure 
based on thin-wire scheme of FDTD to characterize 
transient pulse responses of some communication 
monopole antennas. This paper is organized as follows. 


In section II, the thin-wire model using subcellular 
extension in the orthogonal coordinate system based on 


the method of time-domain finite difference is presented. 


In section II, some numerical investigations are carried 
out to explore the transient effects of an incident EMP 
on some communication antennas, with feed—point 
transient voltage obtained and compared with the 
measurements. Some conclusions are made in Section 
IV. 


Il. ©FROMULATION 


A. The Time-domain finite difference method 


Since the time-domain finite difference (FDTD) is 
proposed by 1966, it developed fast, and has been 
implanted in many fields. The FDTD technique is a 
computationally efficient means of directly solving 
Maxwell’s time domain-dependent curl equations of 
their equivalent integral equations using the finite- 
difference technique. The FDTD discrete the spatial 
domain by the Yee cell and the curl equations 
approximately by the second-order finite difference. 


The EMC/EMI problems can be considered to solve 
the Maxwell equations on the special boundary 
condition, for the opening region problem, the boundary 
condition is the radiation boundary. Time-domain 
Maxwell is as follows: 
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where H is magnetic field distribution, B is magnetic 
flux density distribution, Eis electric field distribution, 
and D is electric flux density distribution. In the 
orthogonal coordinate system, the two curl equations 
can be written as: 
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Based on the Yee cell, we can use a centered finite- 
difference approximation with the second-order 
accuracy with respect to x-, y-, and z-coordinate variable 
or time variable t, and finally obtain the updated 
equations for electric field component EF, , £, and 


H, and H, . For the 


opening region problem, the absorbing boundary 
condition (ABC) is used to truncate the solving domain. 


magnetic field component H 


x? 


B. The equivalent resistive thin-wire FDTD model 


Recently, a resistive thin-wire FDTD model has 
been proposed, and the antenna excited by the coaxial 
lines is simulated by an equivalent model without any 
additional transmission line (TL) cell modeling [10]. 
Now the thin-wire FDTD model is improved to analyze 
the transient response of the linear monopole antenna. 
By the quasi-static approximation, the aperture in Fig.1 
(a) can be represented by the magnetic frill current, and 
according the transmission line theory, the total voltage 
at the aperture equals the superposition of the incident 
and the reflected voltage, respectively, as: 


Vas) =Vine(t) + Vie O 
If the characteristic impedance is Z,, 
voltage at the aperture can be written as: 


(2-4) 
the reflected 


Vie ) = Vine) Zot (2-5) 
Considering the peopacation mode is the TEM mode, 
so the equivalent magnetic frill current can be written as: 


M(x, y,D =[2V,,.O+Z,J@]/[aln(b/a)] (2-6) 


inc 


According to the Ampere’s law, the total current at 
the feed point can be calculated as: 
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Further, the update equation for the magnetic field at 

feed point can be expressed as: 
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The update equation for the magnetic field at the 


other Yee cell in Fig.l (b) can also be obtained 
according the the Ampere’s law as follows: 
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Fig. 1 the adjacent Yee cell and the structure of the monopole 
antenna(a)The equivalent resistive thin-wire FDTD model(b)the 
FDTD updated cell. 
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C. Intentional EMI Exciting Pulse 


According to the experimental date present of the 
high-power wideband electromagnetic pulse source in 
[11], an analytical equation that describes the waveform 
is found as follows. 


6 J 4; B 
nese, (10) 
i=l 
Where the fitting coefficient E, , ¢, and w, are 
understood, but suppressed here. 
On the other hand, we also suppose the incident 
EMP can be a slow-, medium-, and fast-EMP, 
respectively, which described by the two exponential 


functions realized by the parameters % and “ [12]. 
E, (t) =E,k(a, Bye" -e) (11a) 
_pinla)=In(B) —_, Infaz)-In() 
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III. . NUMERICAL RESULTS AND DISCUSSIONS 


Based on the above mathematical, computer code 
are developed for the performing the transient 
electromagnetic responses of the wideband monopole 
antennas illuminated by an external electromagnetic 
pulse with different waveforms and parameters above 
presented. 


Next we considera monopole antenna mounted on 
a finite size ground, shown in Fig. 2. The cylindrical 
antenna is connected to the front-end by the coaxial line. 
The height of the vertical cylinder conductor ish, the 
inner an outer radii of the coaxial line is a andb, 
respectively. The ratio b/a is set to be 2.3, for this 
problem. 


Outer 
Conductor 


Inner 
Conductor 


II 


Coaxial 
Dielectric 


aaa 


Fig. 2. A monopole antenna with coaxial lien illuminated by an 
external EMP. 
The FDTD cell _ size is_ selected as 


Ax = Ay = 1.4605 mm and Ay=0.635 mm. The time 
increment is chosen to satisfy the Courant-Friedrichs- 


Levy (CFL) as At=Ay/2c, , 
C, denotes the light velocity. The unaxial perfect 


matched layer (UPML) is implanted to truncate the 
computation domain [13]. 


condition where 


When the monopole antenna is illuminated by 
external EMP source, the induced current and voltage 
will produce at the feed point of the antenna, and 
further the induced energy will propagate along the 
transmission system into the front-end part, and may 
bring the interfere, even the physical damage to the 
electric system, so it is important to capture the 
transient response of the monopole antenna which 
widely used in shortwave communications systems. 
Figure3 show the induced voltage at the feed point 
when the monopole antenna illuminated by the EMP 
source [10]. From the figure, it is found that the 
magnitude of the induced voltage as high as 19.8 kV. 


LD) 
o 


S The Ultrawide band EMP 
x —e— Induced Voltage 
—~ 10 
o 
io) 
© 
= 
So 
xe) 
o 
is) 
3 -10 
= 

-20 

0 1 2 
Time (ns) 
Fig. 3. The transient response of the induced voltage on the feed 


point. 

Finally, in addition to the high-power wideband 
electromagnetic pulse present in equation (10), we 
further impose an low-, medium-, and fast-EMP of 
E£, =50 kV/m_ illuminated onto the 


monopole antenna, respectively. 


investigated 


Fig.4 shows the transient voltage response on the 
feed point of the monopole antenna in Fig.2, which is 
also illuminated by a high-power fast-, medium-and 
slow-EMP, respectively. It is shown that among three 
types of intentional high-power EMP incidence, the 
fast-EMP excites the strongest coupled voltage in time 
and spectral domain on the feed point of the monopole 
antenna. 
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Fig. 4. The transient response of the induced voltage on the feed 
point. (a)slow-EMP; (b) medium-EMP; and (c) fast-EMP, 
respectively. 


IV. CONCLUSION 


In this paper, a resistive thin-wire FDTD model is 
improve and applied to efficiently simulate the linear 
monopole antenna. The numerical procedure is 
developed for capturing the transient response of 
monopole antennas, and compared with different 
incident source parameters. Various numerical results 
were presented to demonstrate the transient responses 
characteristic of the monopole antenna. This study can 
severe as the first step and input terms for the planning 
of electromagnetic pulse protection of these special 
antennas against the high-power EMP sources. 
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